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Abstract 



We discuss the Donsker-Varadhan theory of large deviations in the frame- 
work of Hamiltonian systems thermostated by a Gaussian stochastic coupling. 
[ We derive a general formula for the Donsker-Varadhan large deviation func- 

tional for dynamics which satisfy natural properties under time reversal. Next, 
■ we discuss the characterization of the stationary states as the solution of a 

, variational principle and its relation to the minimum entropy production prin- 

ciple. Finally, we compute the large deviation functional of the current in the 
case of a harmonic chain thermostated by a Gaussian stochastic coupling. 

(N ' 
> , 

^ 1 Introduction. 

In the recent years, several studies of large systems out of equilibrium through fluc- 
CN ■ tuation theory have been made [21 [Sj HJ [5l [161 HI [HI [2Q]. In a recent series of 

papers [281 O [291 [SOI [M], it has been understood that in random systems driven 
O ■ out of equilibrium, the theory of large deviations provides naturally a variational 

characterization of the steady states which is related to the minimum entropy pro- 
K/i ; duction principle. In this paper, we pursue this approach in the framework of 

' thermostated lattices of Hamiltonian oscillators and investigate in this setting the 

^ Donsker-Varadhan large deviation theory [Til 112] • The existence of a large devia- 

tion principle for a chain of Hamiltonian oscillators has been already established in 
[391 [35] , here we focus on the physical interpretation of the functional in terms of 
entropy production and on exact computations for Gaussian dynamics. 

Thermostated lattices of Hamiltonian oscillators arise naturally in physics as 
they model either the contact with an environment, the nonlinearity of the dynam- 
ics or the randomness of initial conditions by effective stochastic Gaussian terms 
added to Hamilton's equations of motion. For the class of models considered in this 
paper and introduced in section [21 the noise is degenerate as it acts only on some 
coordinates (the momenta). Building on ^24j, we introduce dynamics with a deter- 
ministic driving force which are analogous to the asymmetric dynamics in lattice 
gas models. These models satisfy the generalized detailed balance which is a relation 
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between the generator of the dynamics, a reference measure and a time-reversal 
breaking function that is a hnear combination of the local energy currents. 

In section El we rephrase the Donsker-Varadhan theory in our setting [32| and 
establish a useful formula (1411) for the large deviation functional of processes which 
have natural properties under time reversal. This generalizes the formula expressing 
the functional as the Dirichlet form of the process in the case of reversible dynamics. 

In section 1321 we discuss, along the lines of [211 [291 |30l |31] , the characterization 
of the stationary measures out of equilibrium in terms of a variational principle and 
we stress the differences and the peculiarities inherent to the Hamiltonian struc- 
ture. In particular, Hamiltonian dynamics leaves the Gibbs (or Shannon) entropy 
invariant, thus it comes as no surprise that the so-called minimum (Gibbs) entropy 
production principle fails to determine correctly stationary states in systems out of 
equilibrium or even in equilibrium: it is blind to the Hamiltonian dynamics. The 
only part of the dynamics modifying entropy is the diffusion modeling the action of 
the thermostats which, in our case, act only on the momenta variables. The degen- 
eracy of the dynamics (non-ellipticity) is the main reason for which the minimum 
entropy production fails to select the correct stationary state. We remark that in 
[8j (section 3.3), this principle was shown to be valid in the case of an Hamiltonian 
system coupled to a diffusion acting on all the variables. We analyze the interplay 
between the entropy producing part of the dynamics and the purely Hamiltonian 
part and show that the Donsker-Varadhan functional provides a variational charac- 
terization of the stationary measure (both in equilibrium and out of equilibrium) 
related to the minimum entropy production principle. 

Finally in section 13.41 we recover from the symmetries of the Donsker-Varadhan 
functional the Gallavotti-Cohen symmetry [IHllISlliniEIllSlISSlESllIlES]. 

In section m we focus on the large deviations of the heat current for thermostated 
lattices of harmonic oscillators. We compute exactly the functional and relate it to 
previous expressions derived for lattice gas models [3l HI [5] . 



2 Models 

We will first recall the general framework of Hamiltonian dynamics and define the 
relevant physical quantities in this context. In section [2^ the notion of generalized 
detailed balance is introduced. It will be an important feature of the Hamiltonian 
systems with a stochastic forcing considered in this paper (see sections 12. 3^ 12.41 and 
I23|). 

2.1 Hamiltonian dynamics 

We consider a one-dimensional lattice of N particles {q,p) = {qi,Pi)i<i<N, each one 
moving around an equilibrium position, the position of the i-th particle is denoted 
by Qi and its momentum by pi. The systems we have in mind are described by a 
Hamiltonian, which is the energy function of the set of particles. 



N 



H{q,p) = Y: 



i=l 



Qi) + U{qi- gi_i)) 
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We will consider either periodic boundary conditions (with the convention ^at+i = qi 
and go = Qn) or open systems with the convention go = Qn+i =0. V is the 
potential energy corresponding to an interaction with an external substrate. U 
describes the potential energy of the interaction between nearest-neighbours. Precise 
assumptions on the growth of the potentials will be detailed in Section 13.11 along 
with the mathematical statements. 

The positions qi and momenta pi of the particles obey the Hamilton's equations, 

. dH . dH 

qi = Pi = ■ (2) 

opi dqi 

The generator of the Hamiltonian dynamics is given by 

"^-^ dH d dH d 
r{ dqi dpi dpi dq 



^ dH d dH d 
LH = j:-7^7Zr + lZ-7^- (3) 



In order to describe the propagation of heat in the lattice, one defines a local 
energy function, 

2 

hiip,q) = y + V{qi) + i(f/(g, - g,+i) + f/(g,_i - g,)), (4) 

such that H = J2f=ihi. The local energy transfer is identified with the transfer 
of mechanical energy between nearest-neighbours. The energy current is therefore 
defined through the time evolution of the local energy, 

^ = Lffhi = ji - (5) 
at 

with ji the microscopic current of energy or heat between oscillator i and i + 1 

ji = -\U\qi - qi+i){Pi+Pi+i)- (6) 

In section HI we will investigate the large deviations of the spatial average of the 
current defined by 

1 ^ 

%>p) = ^Ej- (7) 

i=i 

When such a Hamiltonian system is in thermal equilibrium at a temperature 
T = f3~^, its statistical properties are described by the Boltzmann-Gibbs probability 
distribution over the phase space R^^, 

P(P, 2) = ^ exp [-^H{p, q)j . (8) 

Remark that with respect to that distribution, {pf) = T and the averaged cur- 
rent is null (ji) = 0. This last identity is straightforward because the equilibrium 
distribution is even under the reversal of momenta p — > —p while the current is odd. 



2.2 The generalized detailed balance 

We now recall the notion of reversibility in the framework of lattices of coupled 
oscillators. Then, we introduce the notion of generalized detailed balance which will 
apply when the chain of oscillators is forced out of equilibrium by stochastic heat 
baths (sections 12.31 and 1 2.41) or a driving force (section [2751) . 

In Hamiltonian dynamics, the variables p basically describe the velocities of 
the particles and are therefore odd functions under time reversal. When isolated, 
those systems are reversible in the sense that if one lets evolve the particles from 
some initial conditions at time up to some time T and then reverse all momenta, 
Hamilton's dynamics will take back the particles through the same trajectory but 
with reverse momenta. Thus in the case of Hamiltonian dynamics (with possibly a 
stochastic forcing), the standard notion of detailed balance has to take into account 
the reversal of momenta. We denote by H the operator which reverses momenta 

n(g,p) = (g,-p). (9) 

Let Pt (^{q,p), {q',p')^ be the semi-group associated to an Hamiltonian dynamics with 
a stochastic perturbation (for a precise definition see sections 12731 and 12741) with initial 
data {q,p) and final data {q',p') at time t. We consider a probability measure with 
density p{q,p) symmetric wrt to time reversal Hp = p. The dynamics is reversible 
wrt to the density p{q,p) if for any time t 

p {q,p) Pt {iq,p), ii,^)) = p {i,/) Pt (H(g',pO, n(g,p)) . (10) 

Let L be the generator of this dynamics and let L'^ denote the adjoint with respect 
to the Lebesgue measure. The adjoint L* of the operator L with respect to the 
measure with density p{q,p) is defined by 

L;fiq,p) = p-\q,p) L^ip{q,p)f{q,p)) (H) 
for any regular function f{q,p). Alternatively, L* satisfies 

L* is therefore interpreted as the generator of the time reversed dynamics sampled 
with initial data distributed according to p. Applying (ITOl) for an infinitesimal 
amount of time, the equivalent form of the detailed balance relation can be obtained 

l; = ulu, (12) 

where we used that the density p satisfies Hp = p. For stochastic dynamics, the 
usual detailed balance relation does not involve the time reversal operator, however 
we keep the same terminology for simplicity. 

For general Hamiltonian systems coupled to stochastic thermostats, the re- 
versibility may not hold and (fT2i) has to be replaced by the generalized detailed 
balance relation (see for example [H]) which can be defined as follows 
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Definition 1 We consider an evolution with generator L and p{q,p) a reference 
measure such that Up = p. If there exists a function a{q,p) such that 

L* = nLn + a, and Ha = -a , (13) 

then the triplet {L, p, a) is said to satisfy a generalized detailed balance relation. By 
convention in [TS\) . a acts as a multiplication operator. 

This definition will be illustrated in the examples introduced in sections I2.3[ 12.41 
and 12.51 In equilibrium dynamics (when a = 0) and when sampled with initial con- 
ditions distributed according to the equilibrium measure, the dynamics is equivalent 
to its time-reverse, this is the detailed balance principle. When the system is sub- 
ject to some non-equilibrium dynamics forcing a heat current through the system, 
the stationary state looses its invariance under time-reversal (a 7^ 0). In Markovian 
systems, the generalized detailed balance is a central feature of the Gallavotti-Cohen 
symmetry [26l [161 [HI [21 [SSI [Ml [35] for the large deviation functional of the function 
a (see section [3^ . 

Typically a is a function of the microscopic currents of energy in the lattice 
and the reference measure is an equilibrium or a local equilibrium distribution. A 
given dynamics may satisfy a generalized detailed balance relation with respect to 
different reference measures and different functions a (see section 12. 4p . 



2.3 Heat in the bulk 

The simplest perturbation of the Hamiltonian dynamics ([2]) is to couple each oscil- 
lator to a heat bath. The boundary conditions are fixed go = qN+i = and each 
oscillator i = 1, . . . , N evolves according to 

dqi = pidt (14) 
dH 



dPi = -liPidt - -T:—dt + J2-fiTidwi , 
dqi \' 

where the Wi are standard independent Brownian motions, Tj = is the temper- 
ature of each heat bath and 7^ > is the friction. The generator of the dynamics is 
given by 

L = Ls + Lh (15) 
where the Hamiltonian part Lh was introduced in ([3]) and the symmetric part is 

N g q2 

~i opi dpi 

When the temperatures {Tj}j are not equal, the reversibility is lost and the invariant 
measure unknown. The self- consistent chain [6J, for which the temperatures {Tj}j 
are tuned in order to maintain a constant average microscopic current along the 
chain, falls in the framework of the dynamics ([1] 



We check now that this dynamics satisfies the generalized detailed balance (see 
definition [1]). Given a collection of inverse temperatures (3 = {jSi, . . . ,Pn), we take 
as a reference measure the following density 



1 ^ 



P(3{p,q) = i^exp(-H{p,q)] , with H {p, q) = J2 Pi Q) ^ (17) 



i=l 



where hi are the local energies introduced in (jlj) and (3^^ = Ti is the temperature of 
the heat bath at site i. We compute now L*^ the adjoint of L with respect to the 
measure pp. Since Ls = riL^n and the stochastic part of the dynamics is reversible 
with respect to the measure p/?, one has for any function / 

Lsfiq.v) = p^\q,p)Ll{pp{q,p)f{q,p)) . (18) 

We turn now to the Hamiltonian part. Thanks to the relation lJj = —Lh, we get 

N 

pfL\pp_ = Y.[3,LHK. (19) 

i=l 

From ([S]), one has LhK = ji — with the convention Jq = = 0, thus 

N N-1 

pfL\pp_ = E A(J. - 3^-l) = - E - m ■ (20) 
i=i 1=1 

Since HLh^ = —L^ = L^h, this implies that 

Af-l 

L;^=nLn + a^, with a^ = - E(A+i-A)ji- (21) 

i=l 

Thus the triplet {L^pp^ap) satisfies a generalized detailed balance relation and ap 
is a linear combination of the local currents. 



2.4 Heat at the boundary 

If there are no heat baths in the bulk (7j = for i = 2, . . . , A^ — 1), the dynamics ^ 
represents a crystal of atoms heated at two different temperatures at its boundaries. 
The equations of motion are given by, 

dqi = pidt, i = l,...,N, (22) 
dH 



and, 



dpi = - — {p,q)dt, i = 2,...,N -1, 
dqi - - 



OH 



dpi = -——{p,q)dt-^pidt + J2^Tidwi, (23) 
dqi - - " 

'^Pn = — {p,q)dt-'ypNdt+ J2rfTNdwN, 
dqN - - ^ 

where Ti and T/v stand for the temperature of the left and right reservoirs, respec- 
tively, whereas wi and wn are two independent standard Wiener processes. 

When Ti = T/v = T = (3~^, the Gibbs measure (E]) is invariant (stationary) for 
the stochastic dynamics defined above. For two different temperatures, existence, 
uniqueness and exponential convergence to an unique invariant state has been es- 
tablished under fairly general conditions on the potentials U and V p, IT3 | [T5l [M] . 



A computation similar to the case of heat baths in the bulk fl2T|) shows that the 
dynamics fl22|) satisfies the generalized detailed balance (see definition [1]) for any 
reference measure pp of the form f[T7j) 

JV-l 

L;^=nLn + a^, with = - E (/^m - , (24) 

i=l 

provided the collection of inverse temperatures /? = . . . , (3^) is such that /^f^ = 
Ti and I3jj^ = Tjy. Thus any such that /3f ^ = Ti and /^j^"*^ = will satisfy the 
Gallavotti-Cohen symmetry relation (see section [3^ . This was already observed in 

[HEll. 



2.5 Asymmetric periodic chain 

Building on a previous work [23], we introduce new dynamics by adding a mechanical 
force which creates a current through the system. On a periodic lattice, these 
dynamics lead to non-equilibrium systems with non-vanishing currents and in this 
sense, they are reminiscent of the asymmetric processes in lattice gas dynamics [36] . 

We start with an heuristic discussion before giving the definition of the dynamics. 
If the heat baths are at different temperatures, the dynamics fl22l) is no longer 
reversible with respect to the Gibbs measure or any "local equilibrium" measure 
As observed in ( 12^ . the function a which breaks the reversibility is a linear 
combination of the local energy currents. We show that adding an appropriate 
mechanical force allows to modify at will the coefficients of this linear combination. 
In particular, the reversibility may be restored in (!22|) by tuning the intensity of 
the additional mechanical force. From the point of view of the generalized detailed 
balance, the action of the mechanical force is equivalent to the action of the local 
temperature gradient. 

Let us first show that one may choose a (non-Hamiltonian) force which modifies 
the coefficients of the combination of local currents of energy in a. Take as generator 
of the dynamics, 

L = Ls + LH + Le_, (25) 

where Lh is the generator of the Hamiltonian dynamics (I3l), Ls the generator of the 
two stochastic reservoirs (1231) 



Opi Opi OpN Opj^ 

and the contribution of the mechanical drift is given by the antisymmetric operator 



N o 

Le_ = - Y.\W^^lU'{q^^^ - q,) + d,U'{q, - g.+i))^^ , (26) 

i=l '^Pi 



To study the reversibility (lT2l) properties of this dynamics, we first compute, 

N 

P'p^LIpp = -^i(0i_if/'(g,_i-g,) + ^it/'(gi-gi+i))pi 

i=l 



i=l i=l 
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Since the operator Ls + Lh satisfies a generalized detailed balance relation wrt 
= —J2j^Ji^iPi+i ~ Pi)ji (Ell), see that the dynamics fl25l) satisfies now the 
generalized detailed balance (see definition [T]) for the measure p/3 ffTTj) with respect 
to a^^0 

N-l 

L;^ = Um + with ap,e = E (A " A+i + Oi)ji . (27) 

1=1 

When the relation 9i 7^ (A+i ~ A) for some then the system is driven out of 
equilibrium. However, with the choice 6*0 = ^tv = and 9i = (/3j+i— reversibility is 
restored. This is the key point which allows to identify the strength of the mechanical 
force with the action of the local gradient of temperature. 

For 9i = (Pi+i—Pi) and slowly varying temperatures of the form = Ti = T(-^) 
where T is a smooth function, the generator Lg becomes at lowest order in 

= I ^^(^'(*- - *) + - ,.«))|- . (28) 

Note that in the sum over i, the prefactors depending on the temperature and its 
gradient are basically constant when i varies over distances much smaller than N. 

The asymmetric periodic chain [23], described below, is a chain with periodic 
boundary conditions and a dynamics made of three parts. The first one corresponds 
to a usual Langevin dynamics for each oscillator, the second one to the Hamiltonian 
dynamics on the lattice and the third one is the previous generator (!28|) with constant 
prefactors (but arbitrary values). Namely, it is defined as 

L = Ls + Lh + U, (29) 

with, 

i dpi dpi 
Lh is defined in (EJ and the part driving the system out of equilibrium is, 

Lr = mU\q^-l - q^) + U'{q, - • (31) 

In terms of equations of motion the dynamics is described as follows, 

dqi = Pidt (32) 
dH 



dpi = -IPidt - -Q^dt - —{U\qi^i - qi) + U'{qi - qi+i))dt + ^jl^iTdWi , 

where r G R is the new parameter regulating the strength of the non-equilibrium 
forcing and the Wi are standard independent Brownian motion i = 1, . . . , iV. Periodic 
boundary conditions means go = qN and ^at+i = q\. Compared to the dynamics 
(|T^ . the new term proportional to r is the non-equilibrium part of the dynamics. As 
should be clear from the above argument it is responsible for the breaking of the time- 
reversal symmetry of the equilibrium dynamics (r = 0). Indeed, for the dynamics 
), taking the Gibbs measure ([S]) at constant temperature T as the reference 
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measure, the generalized detailed balance holds with a function a proportional to 
the total current ([7]) 

^=j^2^3i = ^J- (33) 

1=1 

We will show in subsection 13 . 31 that when t 7^ 0, t{J)p > for a stationary measure 
p. Thus, the dynamics ensures the existence of an average non-vanishing energy 
current in the stationary state. In section HI the current large deviations for the 
forced dynamics (l32l) with harmonic potentials will be computed. We will see there 
that (132|) are the optimal dynamics to realize current deviations. 



3 The large deviation functional 

The goal of this section is to rephrase the Donsker-Varadhan theory [HI |T2] in the 
framework of the Hamiltonian systems coupled to Gaussian stochastic thermostats 
and to discuss the relation with entropy production. 



3.1 The Donsker-Varadhan functional 

In order to cover all the examples introduced in sections 12.31 and 12.51 we consider 
the general dynamics defined by 

dqi = Pidt (34) 

dPi = -liPidt - —dt - \Ti{ei^iU\qi - Qi^i) + 6iU'{q,+i - qi))dt + J^'^iTidWi , 
dqi ^ ^ 

with 7j > 0,Tj > for % = 1,...,N. We stress the fact that the noise acts at 
each site i. To simplify notation, we restrict to open systems (with the convention 
% = Qn+i = 0, 6*0 = ^at = 0), but similar results hold with periodic boundary 
conditions (with the convention ^at+i = qi and go = Qn)- We denote by P the 
probability of the evolution (!34|) starting from a given initial data (which will play 
no role in the large t asymptotic). 

Typically, we shall be interested in the deviations over time of some physical 
quantities like the heat current ([7]). Let the empirical distribution ut be defined by, 

Ut{A) = ] f\A{{p{s),q{s)))ds, 



t Jo 

where 1^ is the indicator function of a set A C R^^. If the dynamics is ergodic 
i^t converges to the ergodic invariant measure. We now look at the asymptotic 
probability P[z/f ~ /i] that the empirical distribution is close to the distribution of 
a given measure /z (in the sense of the weak convergence topology [39]) for large t. 
Under suitable hypothesis (see Proposition [1]) , the dynamics obeys a large deviation 
principle with rate function / and asymptotically in t 

P[z/i~/i]~exp(-t/(/i)). (35) 

Furthermore, from the Donsker-Varadhan theory [Til [l2l [39] the functional is given 

by 

/(/x) = sup I - /^^\ ge Cr(R'^; [1, oo[) I , (36) 
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9 /. 



where C^(R^^; [l,C)o[) is the set of bounded infinitely differentiable functions in 
R^^ taking values larger or equal to 1. We refer to Lemma 6.3.7 of [L0\ for the 
variational expression of / in the case of hypoelliptic diffusions. 

The following Proposition justifies the validity of the large deviation principle 

Proposition 1 Suppose that the potentials V and U of the Hamiltonian (QP are 
convex, twice differentiable and satisfy 

V"iq)>5 and E Viq.) + ^^^^ " ^^"^^ t ^^^^^^ " > ^ E {U'iq. - , 

(37) 

for some constant 5 > 0. 

The dynamics (34\ ) (with Wi,'yi > 0) obeys a large deviation principle with a 
functional I given by the Donsker-Varadhan theory provided max j 16*41 < Tq 
and maxj |Tj — Tj+il < Aq, where Aq and tq are two constants depending only on 
S, {ji, Tj}j. Furthermore, the previous assumptions on the potentials ensure that 
the current ^ is exponentially integrable: for any A small enough 

lim sup - log P ^exp ^A J dsJ{p{s),q{s))^^ < oo . (38) 

The proof heavily relies on previous results in the paper [39] and it is postponed 
to the Appendix. The assumptions of Proposition [1] on the potentials V and U are 
not optimal, but they are sufficient to cover a wide class of physical examples. In 
particular, similar statements hold also for any local modifications of the potentials 
V and U. Remark that when the reservoirs act only at the boundary (7j = 0, z G 
{2, . . . ,N — 1}), the large deviation principle for the current has been justified in 
|35] (for a different class of Hamiltonians) . 

We are going to rewrite the functional / (!36|) in a more explicit form. The 
generator of the dynamics can be decomposed as L = Ls+L^, with a symmetric 
part Ls (p^Gjl due to the noise and an antisymmetric part La = h{L — IILn). We 
introduce the notation, 

N 

T{f,g) = 2j2^m{dpJ){dp^g), (39) 

i=l 

for any smooth functions /, g. 

Proposition 2 Let p be a measure and c^q^p) a function such that 

l; = nLn + a . (4o) 

Let p be a measure absolutely continuous with respect to the measure p, with f = 
dp/dp G C°°(R2^;R) such that (HVpT/lh) < oo. Then the functional I ^ is 

given by 

Hp) = - (f^Lsf-) + Kip) - i (a)^ , (41) 



in 



with 

Kifi) = - inf Q {TiW, W))^ + \ {LaW)}j > , (42) 

where the infimum is taken over the smooth functions W G C°°(R^^; R) such that 
W e L2(/i) and \VpW\ G L^ifi). 

From (^71) . we see that the dynamics satisfies the generahzed detailed bal- 
ance for the measure p/^ f|T71) . the assumption pUj) of Proposition [2] is satisfied with 
and a^^e 

L*^ = IILn + ap^e, with cr^^g = ~ A+i + • 

i=l 

Remark that the generalized detailed balance requires symmetry assumptions on 
the reference measure and on a which are not necessary for the Proposition [2] to 
hold. 

Proof. Let / = dp/ dp and g E C^iR"^^; [1, oo[), then, 

[g-'Lg) = (fg-'Lg) . (43) 



Write now in (H3!) . g = \fhf for some h > and symmetrize the expression with 
the help of the adjoint operator of L, 

(fg-'Lg)^ = \ (^(^fih^iLihif^))^ + (^fih^L*ih~ifi))^^ (44) 

We note the relation, 

L{(Pij) =(PLij + ijL^ + T{i:,^), (45) 
And analogously for L* = IILn + a = Ls — La + cr (by hypothesis), 

L;{(f)^) = <p{Ls - La)^ + ^{Ls - La)<P + r{^, 0) + CT<I)^. (46) 

Applying those formulas for L (resp. L*) with (p = and ip = (resp. (f) = h'"^ 

and ip = f^) using the generalized detailed balance relation L* = ULU + a and 
computing systematically all derivatives, one gets. 



P 8 \ /i2 / . - \ ' 



(47) 



Writing = e^ and using the definition (136]) of /(/x), one finally gets the variational 
formula restricted to functions W such that {W + log/) G C^{R'^^;R). To 
extend the class of functions W in the variational formula (H2l) . we consider an 
appropriate sequence of functions ipn ^ C°°(R^^;R) such that ■?/'„ = 1 in the ball 
of radius n and = outside the ball of radius n + 1. Then the sequence Wn = 
Wipn - (1 - '0n)log/ approximates any function W G C°°(R^^;R) fl L'^{p) with 
|Viy| G L2(/i). 

The fact that K{p) > follows by choosing = in fH2l). 
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One of the important features of the representation fj4T|) for /(/x) is the presence 
of K{fi). As we shall see below, in Proposition [3] and in the example (l52l) . it is 
in general neither infinite or zero even in the case of dynamics satisfying detailed 
balance, i.e for equilibrium dynamics. K{n) gathers the thermalizing effect of the 
noise which is transmitted from the p variables to q variables. It corresponds to the 
"traffic" in the terminology of ^29j. It is invariant under reversal of the sign of the 
momenta (i.e K{fi) = K{Iln)) and when non-equilibrium forces are included it is 
invariant under a change of sign of the non-equilibrium parameter as we shall see 
in the gaussian systems of the last section. It measures the dynamical "activity" in 
the Hamiltonian system irrespective of the sign of the non-equilibrium parameter. 
We shall see its role in the variational characterization of stationary states in section 
13.21 We now compute a more explicit form for K[ii). 

Proposition 3 Let ^ be a measure on R^^ with smooth density wrt the Lebesgue 
measure dfi{p,q) = exp{—^)dpdq. IfW& -^^(/^) is a solution of 

N 

Y: -l^T^^lW + 7^T. d,^^d,^W = -La^ , (48) 

with \VpW\ e L'^ifi) then 

K{^^) = ^ {t{W,W))^ = -\ (LaW)^ . (49) 

Moreover, (t{W,W)^ is independent from the solution of ( f^g] ). Finally, K is 
symmetric wrt time reversal K{J1^) = K{fi). 

Proof. Starting from the definition (142!) of K{fi) and integrating by parts, leads to 
^ (r(l^, W))^ = ^Y: -l^T^ {Wdlw)^ + {Wd,^^d,W)^ . (50) 

i 

Integrating also by parts the second term in P2l) gives, 

\{LaW)^ = \{WLa^)^. (51) 

Therefore, the solution to the variational problem in (142!) gives the equation (HHj) . 
Using then ([50D and ^ and the definition (gSD of K{ii) finally yields (09]). The 
second equality in (H^ follows by combining (150]) and (HSj) . To see that K{jj) does 
not depend on the solution to (H5!) . let W be two different solutions, then, 

N 

Y: -l^Ti^l{w -w) + ^,w,^^d,xw - ly) = 0. 

i=l 

Multiplying this by -|- and integrating by parts with respect to /i, yields 

(V{W + W,W -W)) ^ = Q, 

and thus, {V{W, W)) ^ = (t{W, W)) . 
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The measure Ufi has a density given by exp(— 11$) and — nVT is the correspond- 
ing solution of (HHj) . Thus, the symmetry of K follows from the identity 

This concludes the proof. 



Observe also that if is such that L^^ = then K{^) = 0. We want to 
emphasize that in general K is non-zero and we will illustrate this in the simple case 
of a one dimensional harmonic oscillator coupled to a heat bath at temperature T 

d d d 

L = -IP-^ + iT— - uj^qTT + Pt- , (52) 
op op"^ op oq 

with Ls = —lP-§^ + iT-^ and La = —co'^q^ ~^ P'§^- ^^^^ compute K{jj,), where 
II is the Gaussian measure 

/^(?>P) = ,^ ^ A,^ ,1 exp (-{a^ + buj'^'^)) . 
{2Txhuj^Y{2Tia)^ V ^ ^ / 



In that case, (148|) becomes, 

-^TdlW + ^TapdpW = u'^pqib - a) . 
A solution is readily found, 

T77 b - a 2 
W = — —u pq , 
a'yl 

and then, from Proposition [3l 

From this example, we see that K[ijl) vanishes only for Gaussian measures in which 
the average potential energy and kinetic energy are equal, i.e when equipartition of 
energy is realized. In a sense, K measures how much the Gaussian noise acting on 
the momenta variables p is transmitted to the positions variables q. 



3.2 Entropy production 

In this section, we show that the minimum entropy production principle does not 
apply for Hamiltonian dynamics and that the large deviation functional fHTj) is a 
natural extension for a variational characterization of the steady state. Connections 
between large deviation functionals in stochastic systems and entropy production 
were initiated in [28] . In order to identify the average entropy production of a 
dynamics in a given measure, we proceed as in [36] in the context of interacting 
particle systems. Namely, we define the entropy production as the difference between 
the variation of the Gibbs (Shannon) entropy and the transfer of heat by unit time 
due to the action of the thermostats. A similar computation was performed in [27] for 



heat conduction networks and the identification of the average entropy production 
with the Dirichlet form of the process was obtained there and in |14j . 

We consider the dynamics flM|) for which the generator is given by L = L5 + La 
with a symmetric part given by 

dpi 

We first introduce the entropy production. In thermodynamics, the entropy varia- 
tion rate, or entropy production is the transfer of energy per unit time divided by 
the temperature at which the transfer takes place. Therefore, it is natural to define 
the average (with respect to a given measure /i) entropy fiux associated with the 
exchange of energy with the external heat baths as 



where 

i.-g..r.-V.^ + g.,^. (54) 

We have simply divided the contribution of each bath by its temperature Tj: when 
the temperatures are all equal to T then L5 = ^Ls- Computation yields, 

N 2 

L5i^ = E7.(l-|), 
i=i -'^ 

and thus the average entropy fiux due to the coupling to the heat baths in the 
measure fi is. 



N 



^e.t(/^) = E7.((l-^)) • (55) 



i=l 



The Gibbs (or Shannon) entropy for any measure /i is, 

S{ii) = - Jdxflogf. (56) 



We define the entropy production in the chain as 



s(/^) = ^S{^it)\i,t=t. - (^ext{lj) ■ (57) 



We show now that with this definition, the entropy production is always positive. 
Note that, 

d 



— 5'(/ij) = - j dxftLlogft. 



As L = Ls + La and La is a first-order differential operator such that = —La, 
we see that 



d_ 
di 



SilJ-t) = - J dxft{Ls + La) log/* = - J dxftLs log/i. 
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Comparing the density ft with the reference measure pp (see (1171) ) 



= ^ exp ^- J2 A^i^ , (58) 
where Pi = T^^^ for every i such that 7i 7^ (other Pi are arbitrary), we get, 
:S{fit) = -J dxftLslogj^ + J dxftLs (^^Pih^j 

£ N 2 

= - [ dxftLshg^+Y.^, [ ftil-^). (59) 
PfS i=i J J-i 



d 

di' 



We recall that r(/, g) = 2 YliLi liTi dp^/dp^g- For the first term, we note the identity. 



1 



Lsilogh) = h~'Lsh - \h-'r{h, h). 

Combining (155|) . (l59l) and using the fact that (LsYpfs = (because Pi = Tf^ for 
every i such that 7i 7^ 0), with ([57]), this yields 

s(/i) = I / rfxp^(^)r f ^, A=2 [ dxppT I , I , (60) 

with / the density of the measure p. Since r(/, /) > 0, it is easy to see that s{p) > 
and that the infimum is reached when / = p^. As already observed in [27], it is 
also obvious that, as F only involves derivatives with respect to the variables p, one 
can add any function of q in the exponential defining pp. Therefore, even in the 
case of an equilibrium dynamics, when Tj = T, Vi, such that 7, 7^ in L5, the 
minimum entropy production principle does not single out the equilibrium measure 
pj3. This comes from fl58]) which expresses the fact that the Gibbs entropy is invariant 
under the Hamiltonian evolution. Entropy is produced solely by the action of the 
thermostats which act only on the p variables. 

As observed in [28] , the large deviation functional / provides a natural variational 
characterization of the stationary measure, which as we will see below, generalizes 
the minimum production entropy principle. As far as the variational principle is 
concerned, the key observation is that I{p) > and I{p) = if and only if p 
is a stationary measure for the process associated to the generator L (Theorem 
4.2.39 of [ID])- We apply now Proposition [2] to the dynamics (1341) which satisfies 
a generalized detailed balance (I27|) wrt the reference measure p = pp (1581) and 

(yp_,e = Ef=i\Pi - P^+l + 0i)3^. Thus dH]) reads 

I{p) = -^s{p) + K{p) - \ (or^,,)^ , (61) 

where we have used the following identity obtained by integration by parts from 

1 



.s{p) = -l^f^Lsf^) , (62) 

4 \ ' P/3 



for the reference measure pj^. The first term in fl6Tl) is identified with the entropy 
production s of the measure /x, and as we have seen in the previous section, the sec- 
ond one K records the couphng between the positions and momenta. The last term 
in fl^Tl) is the time-reversal breaking term given by combination of the microscopic 
currents. For equilibrium dynamics, by definition, the term a is absent, and the 
presence of K ensures that the minimization of the sum of the first two terms in 
(!6Tll yields the stationary measure univocally. Indeed, as we have explained above, 
minimizing entropy production alone is not sufficient to determine the equilibrium 
distribution. The example (l52l) of a single harmonic oscillator coupled to a heat 
bath is interesting in that respect and we come back briefly to it. We have already 
computed K{ii) in fl53l) and compute now the entropy production We consider 
the reference measure, 

P(P, 9) = ^ exp (^-^(y + ^^|-)) ' (63) 
and a Gaussian measure 

(27r6u;2)^(27ra)^ V ^ ^ / 

Then s{p) = a~^'-)T{a — (3)'^. Combining this with the expression of K{p) ( l53l) . we 

In this simple example, we see that the role of K is to ensure the equipartion of 
energy by the action of the Hamiltonian dynamics, the role of the entropy production 
is to fix the temperature of the equilibrium distribution. 

3.3 Positivity of the energy current in the asymmetric chain. 

We come back to the asymmetric chain defined in fl32l) and use the previous results 
on the representation of the Donsker-Varadhan functional to prove that when the 
parameter r is different from zero then in the stationary state p, r(J)p > 0. We 
note that the hypoellipticity of the process defined by the asymmetric chain may 
be shown by checking the Hormander condition on the generator as in [3lj. For the 
asymmetric periodic chain, no coupling between nearest-neighbour is even required 
because the noise acts on every particle. Hypoellipticity of the process implies the 
smoothness of the probability transition and thus that the stationary state, whenever 
it exists, is described by a smooth density. Irreducibility properties of the process 
may be checked by a control argument and Stroock-Varadhan support theorem [38] 
as explained in [9l [31] . Hypoellipticity and irreducibility imply together that there 
is at most one stationary measure. The existence of a (unique) stationary measure 
follows from those two properties and the existence of a Lyapunov function as in 
( I102p in the appendix (see for instance Theorem 8.7 of [33j). 

Proposition 4 Let U and V be potentials satisfying the hypothesis of proposition^ 
and such that U" is not uniformly equal to 0. Let p be the unique stationary state 
for the dynamics l[2^) . then t{J)p > 0. //r 7^ 0, then t{J)p > 0. 



Proof. Take as a reference measure p/3 = exp(— it was noted in fl33|) that 
the dynamics satisfies a generahzed detailed balance relation with respect to pp and 
a, with a = ^J, T = Write the formula fHT]) of Proposition [2] under the form 

(JSH), 

/(//) = ^s(/i) + i^(p)-i(a)^, (65) 

for a smooth measure /x. As p is a smooth stationary measure, then I{p) = 0. Since 
s{p) and K{p) are positive, one concludes that r(J)^ > 0, using fl65|) . To show that 
if r 7^ 0, (J)p > 0, we proceed by contradiction. Assume then that r 7^ and 
{J)p = 0, then simultaneously s{p) = and K{p) = 0, since the two are always 
non-negative. Writing p = pi3exp{—(p), the first condition implies by construction 
of s(p), (see deOD) that, 

dp^ifi = 0,yi. (66) 

On the other hand, K{p) = implies that the solution W in fH9l) is such that 
dp^W = , Vi. Therefore, using ( HHI) . (y9 must be such that, 

LA{f3H + y^)=0 (67) 

Since La = Lh + and combining ( l25l) with ( 166|) . satisfies 

i=i ^-^ i=i 

with \E'j = f/'(gj_i — gj) + f/'(Q'« ~ Q'i+i)- As this equation holds for any vector p and 
because by (!66l) . does not depend on it means that 

^..^ = . Vz. (69) 

But since (9g-\E'j+i = f/"(Q'i ^ li+i) ^"^^ dq^^-^^i = —U"{qi — qi+i), there cannot be any 
solution to fl69|) . unless U"{q) = for all q. Therefore, we conclude that (J)^ 7^ 
and thus t{J)p > 0, because r(J)p > 0. 

3.4 The Gallavotti-Cohen symmetry 

An interesting aspect of (14T!) in Proposition [21 

/(p) = -(/iL5/n +K{p)-Ua)^, (70) 



is that when the triple {L, p, a) satisfies a generalized detailed balance relation, the 
three terms are either odd or even under the reversal of momenta 11. The fact that 
K{Ilp) = K{p) was established in Proposition [3l In the previous section, we have 
identified the first term in (1701) to the entropy production s{p). It is also invari- 
ant under time-reversal (i.e s(Jlp) = s{p)) when Up = p because by construction 
IIL^n = Ls- And, since a is odd under time reversal (Her = — cr), one may therefore 
write, 

/(n/i) = - (f^Lsf-)^ + Kip) + i ia)^ = lip) + ia)^ . (71) 

Heuristically, (I7T1) implies the Gallavotti-Cohen symmetry for the large deviations 
of the function a integrated over time 



lim — log P 

t^oo t 



1 /■* 

- / (T(s)ds = cr 
t Jo ^ ^ 
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with 

/(a)=inf{/(/x)|(a)^ = a}. (72) 

Since Her = —a, 

ii-a) = inf{/(/x)|(a)^ = -a} = inf{/(^)|(na)^ = a} 
= inf{/(/x)| (a)n^ = a} = inf{/(n/i)| (a)^ = a} . 

Thus one gets from (ITTjl 

/(a) - /(-a) = -a. (73) 

This is the celebrated Gallavotti-Cohen symmetry relation [T8l [T9l [T6l[2T| [221I231I26] . 
For a chain of oscillators in contact with two heat baths, the Gallavotti-Cohen 
symmetry has been rigorously justified in [HI [35]. 



4 Large deviations of the current 

Using the Donsker-Varadhan variational principle (HTl) . we are going derive explicit 
expressions for the large deviation function in the case of harmonic interaction when 
the system size diverges. We consider the asymmetric dynamics fl32|) on a periodic 
chain of length where each oscillator is coupled to a heat bath at temperature T 

dqi = pidt (74) 
dPi = -IPidt - -Q^dt + — (gi+i - qi-i)dt + '^2-fTdwi, 
with harmonic potential 

TV 2 2 2 

H{g,p)=^^J + ^q^ + Y(9^-q.^lY, (75) 
and a pinning ujq > 0. 

We will study the large deviations of the spatial average of the total current ([7]) 

^2 N 

J^l^P) = -^^i^i+i - Qi){Pi+Pi+i)- (76) 

i=l 

The large deviation principle (Proposition [1]) implies that for A and r small enough 
J^Ij, = ^hm ^logP (exp (iVA |J ds J(g(s),p(s)))) , (77) 

where 

Tl^ = snp{J^lM} with J^IM = NX{J)-Fifi), (78) 



where the upper-script r emphasizes the dependency on the asymmetry. We recall 
from Proposition [2] that /^(/i) can be written, 

r(^) = ^.(^) + ir-(/i)-^(j)^, (79) 

where s(/i) denotes the entropy production (l62ll . 

For the harmonic chain (17^ . Gaussian computations lead to the following asymp- 
totics for ^Xjv 
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Proposition 5 

(i) For T = and fixed N, one has for A small 



N rp2, A / L \ 2 rp4, ,8 rp4 8 / i, \ 4 



^0) 



wit/i u;^ = cug + 4cj^ sin (27r;|:)^. 

(^zzj Fzx uo > and two parameters X',t' G R. For large N, with the scaling 
A = X' /N and r = t' /N, one gets 



lim NTl^ = KiXV + iX'YT') , 



where k is the conductivity of the model 

1 f\ u;^sin2(27rx) 

K = - , n/ — zdx. 82 

tug + 4cj2sin2(27rx) ^ ' 

Before deriving Proposition [5l we first comment on the results. Taking the 
Legendre transform of flHTj) shows that for a weak drift r = ^, the large deviations 
of the currrent are given by 

where /at scales, for large A^, like 

In the stationary regime, the mean current has been already computed in terms of 
the conductivity ([82]) (see [21]) 



lim lim { J {q{t) , p{t))) = kt' . 

N—>co t^oo \ — — / 

The expression (!83]l is new in the framework of thermostated Hamiltonian systems: 
it relates the linear response theory [2l| to the large deviations and shows that the 
atypical behavior can still be described by the conductivity k. The reason is that the 
local equilibrium is preserved for small shifts of the current (remark that the fourth 
order term in flHOl) corresponds to corrections to local equilibrium). For weak current 
deviations, fl83]) has a Gaussian structure similar to the one obtained in lattice gases 
with a weak asymmetry [5] [3]. The variance in ( l83l) (given by 2kT^) is related to 
the response coefficient k according to the usual Einstein fluctuation-dissipation. In 
the derivation of ( IHOi) . we will see that the corrections due to K'^^^ arise only at 
the order A^. The functional is only quadratic asymptotically in A^. The Gaussian 
structure will not remain for non vanishing drifts r or larger current deviations. A 
similar behavior can be observed in lattice gases like the weakly/totally asymmetric 
exclusion processes [5]. 

Proof of Proposition O 

1Q 



As the system is periodic and the functional J^x^q is convex, one can look for 
a maximizer of the variational principle fl78l) among the translation-invariant mea- 
sures. The current J in ([7]) is quadratic and the evolution equations fITIl) linear, thus 
the minimum in (1751) will be reached for Gaussian measures (to see this, one can 
check that the maximizer of (178|) is obtained in terms of the principle eigenvector of 
the operator L — \J and the principle eigenvector of its adjoint). Thus it remains to 
optimize the variational principle (|78ll over Gaussian translation-invariant measures. 
We state first Gaussian estimates before deriving (i) and (ii). 

Gaussian estimates: 

For periodic systems, the Gaussian computations are easier in Fourier coordi- 
nates. We set 

„ i X ^ ,-27ri,- „ i 



For notational simplicity, we choose an odd length = 2?2 + 1 and the Fourier 
modes k G {— n, n}. This choice leads to symmetric relations: Pk = P-k and thus, 
|Pfcp = P^kPk and similarly for the variables Qk- 

We will first provide explicit expressions of /^(/i) for a general Gaussian measure 
/i. Given the reference measure po 

Po = e^v{^-^Y.\Pk? + u:l\Qk\^ ' ^^^^ 

the most general form for a Gaussian translation-invariant measure in the variables 
{P.,Q) may be written p = exp(— |;\E')po, with 

n 

^ = E (i (^kPkQ-k + bk[\Pk\^ + col \Qk\^] + Ck[\Pk\^ - col \Qk\^]) (86) 

k=—n 

with Gk = 3?(afc) + iQ'(afc) = —a-k G C and bk = b^k G R, = c^k G R- The 
expectation under p will be denoted < ■ >^, and for each mode k G {—n, n} 

Pk\^) =Skil + 2{bk-Ck))ujlT, {\Qk\^) = Sk{l + 2{bk + Ck))T (87) 
i {PkQ-k - P-kQk)^ = -2^{ak)6kT, {PkQ-k + P-kQk) ^ = -2'<s{ak)5kT 

with = (1 + 2bk)'^ujl - Aclul - \ak\^. 

Using now Proposition [3], we are going to compute K'^{fi) for the dynamics 
and fi as above. In Fourier coordinates the antisymmetric part of the generator 
La = Lh + Lr (??) is given by 



with ujI = uJq+ sin i^ixj^ and 



i. = -.^E™>(^)Q.g^. (89) 



90 



We first show that if Q{ak) 7^ for some k, then K'^{n) = +00. Remark that 
LAlQkl"^ = PkQ~k + P-kQk and that according to fl87|) . (L^IQfcP)^ = if and only if 
^{ttk) = 0. From the definition of K'^ifJ^) fH21) . if $5(0^) 7^ for some k, it is possible 
to choose a sequence of real numbers Cn such that for Wn = CnlQfcP 

Jim (1 (r(iy„, iy„))^ + i {LaW^)^ = \im (P.g-. + P-kQk), = -00 , 



and thus K'^^jj) = +00. So from now on, we assume that is real and that 
afc = —a_k. In that case a solution of (HSj) is 



(90) 

By Proposition [3l one gets, 

1^4(1. 2(.. + c.),(^. n (24)% . (91) 

with ujI = uI + Aoj'^ sin (27r^)^ and 5^^^ = (1 + 2bkfujl - Aclul - a\. 

Note that K'^i^^x) is made of two parts: an equilibrium one (independent of the 
forcing r) and of a non-equilibrium one. 

The entropy production ( 160|) is given by 
^ « / 9vl> 9vi> \ 

= ^ ^ 2iafc(6fc + Cfc) (Pfcg_fe - P-kQk), + «fc (iQfel')^ + 4(&fc + c,)2 (1 )^ 
For Ofc real, one gets from (IHTl) 

n 

s(/i)=7 E 5fc[^^(^~2(&A. + Cfc))+4(6fc + Cfe)2a;2(l + 2(6fc-c,))] . (92) 

fc=— n 

Finally, the current can be represented in Fourier coordinates 
= -i E sin [^^^ {PkQ-k - P-kQk) , 

k=—n \ / 

and from ( 1871) . its expectation is 

uj'^T " ( k\ 

k=—n \ / 



Proof of (i) : 

The functional (l78l) reads for r = 



^,V(^) = -^.(/x) - K%^i) + NX ( J)^ 
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Restricting to translation invariant Gaussian measures fl86|) . the previous Gaussian 
computation lead to 



k=—n 

7 



XTuj sin 211- 



N 



7(1 + 2(6fc + Cfc)) 



(94) 



^ (al{l - 2{hk + Cfc)) + A{hk + Ckfujl{l + 2(6fc - c^))) 
Optimizing over each Fourier mode leads at the order 



Ck 



22V 

7 



sin (27r4 I A + 0(A3 



sin 27r- 



+ 



sin (2. A) A^ 



(95) 



The previous computations imply that at the order A the solution of the variational 
principle (j78|) is given by 



1 



exp 



'2AAr 

. 7 



J{(i,p) Po{q,p) ■ 



(96) 



Finally, asymptotics (ISO]) follow from (195]) . 



Proof of (ii) : 

We observe that flTHj) may be rewritten 



sup^,V(/i) = sup [-r+^^^^^) - K^(/.) + K^^''^\ii)] . (97) 

The point of writing the variational problem under this form is that minimizing only 
JT+2AT easy. The minimizer is the stationary measure of the asymmetric periodic 
chain with non-equilibrium parameter r replaced by r + 2AT^. It was derived in 



Pr+2xMq,p) = ^exp[ --H{q,p) + -^(r + 2\T')J{q,p) . (98) 



For small A and r, we will see that Pr+2XT^ is a good approximation of the minimizer 
of the variational principle J-'lj^ at the lowest order (see also (1961) ). 

jT+2\T'^ is a convex function of the parameters {ak,bk,Ck)- Since Uq > 0, the 
contributions of each modes are independent with an hessian uniformly bounded 
from below in a neighborhood of the minimizer Pt+2\t^- Any measure of the form 
/i = exp(— |;\E')pt-_|_2at2, with as in fl86l) . has a contribution 

r-'''^\p)>C± {al + bl + d} (99) 

k=—n 

for some constant C > independent of N and (a^, 6^5 c^) close to 0. 



9.9. 



From the exact expression fl9T|) . we first evaluate the difference in fl97|) for a 
general Gaussian measure (|86|) 



- K^i^) = ((r + 2AT2)2 - r^)! E 4(1 + 2(6. + c.))|^ sin (^27r^J 

( h\'^ 

(100) 



(Ar + A^T^)! ^5,(1 + 2(6, + Cfc))u;Sin (^27r^^ . 



For large A^, with the scaling A = A'/A^ and r = t' /N, the contribution of the 
/jfi_Pourier mode is at most of order sin (^27r-|:j . From fl^^ . the corrections 
{ak,bk,Ck) have to remain much smaller than ;^sin (2njj^, otherwise the contribu- 
tion of I'^+'^^T^ would be too large as compared to K'^^^^'^^ — K'^. This shows that 
at the lowest order Pt+2\t^ approximates correctly the minimizer of J^x,n- Plugging 
the expression ( l98l) of Pt+2XT^ in (EZD and then in f llOOp gives at the lowest order in 
1/N, 

, X (AV + (A')'T2) CU^ „ 1 ( kV (\ 

For large A^, the last sum may be approximated by an integral 

•^aV(p.+2atO = ^(AV + {XTt') + o (^) , (101) 
where k is the conductivity of the model 



Remark. Note that a direct computation of the largest eigenvalue of the operator 
L — \ J would have led also to J^X.Af- 

5 Conclusion 

We have analyzed the large deviations for thermostated lattices of Hamiltonian 
oscillators by using the Donsker-Varadhan large deviation theory [Til [12] • After 
deriving the large deviation principle for a large class of models (Proposition [1]), we 
analyzed the Donsker-Varadhan functional (HTj) and showed in section 13.21 that it 
provides a variational characterization of the stationary measure which generalizes 
the minimum entropy production principle. 

The current large deviation function has been computed exactly (Proposition 
|5]) for thermostated lattices of harmonic oscillators on a ring. The central issue in 
this computation is the evaluation of the "dynamical activity" K of the station- 
ary measures of systems out of equilibrium. In this paper, we used the fact that 
the explicit form of this measure is known in the model treated. It would be in- 
teresting to generalize these computations to Gaussian self-consistent reservoirs for 
which the temperatures of the heat baths vary spatially. Except for the lack of 
translation invariance, this model is very similar to the one considered here. A very 
challenging issue would be to investigate the current large deviations of the chain of 
non-harmonic oscillators. Our analysis of the Donsker-Varadhan functional applies 



to this type of dynamics, but to derive quantitative estimates of the large deviations 
of the current, some type of approximation seems necessary. In the approximation 
[211 [251 [D El EZ] , the action of the nonhnearity of the dynamics is essentially reduced 
to the action of Gaussian stochastic heat baths, thus one should be able to proceed 
in similar fashion as in the model studied in the previous section. 



Appendix 

In this appendix, we prove Proposition [1] for the dynamics (IMl) . The generator of 
this dynamics is given by L = L5 + + Lq_, with 

fr[ dpi dpi fr[ dqi dpi dpi dqi 

N Q 

Le_ = -Y.\T^{0^-lU'{q^-l-qi) + e,U\q^-q^+l)) — , with 6^ = 6^ " 



and open boundary conditions (go = ^a^+i = 0). The case of periodic boundary 
conditions can be treated similarly. 

The proof heavily relies on corollary 2.2 of the paper j39] where it is proven that 
the large deviation principle applies provided that one can find a function \Ef > 1 
such that the function 

*fe^) = -^(E'^) (102) 

diverges at infinity. This can be understood as follows. From (11021) . we deduce that 
^> is an eigenvector for the operator L + $, thus the Feynman-Kac formula implies 
that for any time t and initial data (p(0), g(0)) 

^(p(0),g(0)) = P (^^(p(t),g(t))exp ds^{p{s),q{s)) 

> P (^exp (^^* ds^p{s),q{s))Yj (103) 

where we used that \E' > 1 in the last inequality. This ensures the exponential 
integrability of /g ds^{p{s),q{s)) uniformly in time. Since $ diverges at infinity, 
we recover tightness properties and the existence of an invariant measure for the 
dynamics (see Theorem 3.1 in [39]). 

Proof. Following [39], we introduce \E'(p, q) = exp{F(p, q) — inf F) with 

1 ^ ^ 
F=-H{p,q) + hY,q,Pi, (104) 

with h a small constant and H given by (see ([1] 



H{p,q) = 2^^ + Y.^iqi) + ^ 



9A 



From fll02p . we get 



N 



^p,q) = -LF{p,q) - J] 7^^.(5,^ F) 2. (105) 



i=l 

We compute each terms 

N N-l 



7' / 1 1 \ 

i=l i=l 

^ 1 

1=1 

where LhH was already computed in (120!) . One has also 

T- 

^^iPi = -liPiqi+p1 - qidg^H{q) - -j{9iU'{qi - g^+i) + 9i-iU'{qi^i - qi))qi 
{dnFf=[-^Pr + hq?j . 



Combining the previous results leads to 

= E (4I - fe) + hq,d,Mq) - b\T,q^ - | 

+ l{^^U'{q^ - q^+l) + O.^^U' {q,^i - qi)) (hT^qi + | 

1 1 1 \ 

+ 4E UT-^p'(9m-9.)fe+Pm), (106) 

Since U and \^ are convex, one gets 

E g.a,/7(g) > E I/(g.) + ' ^^^^^ + ^^^^^^ ' - V{^) - U{0) . 
i=i i=i ^ 

Combining this inequality and fll06p . we get 

Hp,q) > E f ^ - P' + hV{qi) + hU{q, - g,+i) - 6V(0) - 6f/(0) - | 

A 



-62(7.r+ + r+2)g,2 - (^^ + (f/'(g.+i - g.)' +pD , (107) 

where A = maxj |Tj — Tj_|_i|, 7^ = maxj Tj, 71 = maxj Tj and 9 = maxj 

From the assumptions (P7|) on the potentials and fll07p . we see that as soon 

then 

TV 

HP,q) > E + Viq,) + U{q, - g,+i)] - , (109) 
j=i 



for some constants c\ > 0, C2. This is enough to conclude that $ diverges at infinity 
and that the large deviation principle follows from fllOSp (see |39]). 



From fll09p . one can also derive an upper bound on the total current ([7j) 

^p,q) + C2>aN\j{q,p)\, (110) 

for some constant a > 0. From fllU3p . we deduce the exponential bound fl55]) on the 
current. 
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